library(sasld)

# -------------------------------------------------------------------
# ATTENZIONE: la funzione sample e' stata modificata 
# a partire da R 3.6.0 per cui i "vecchi risultati"
# non posso essere riprodotti di default
# http://hostingwin.unitn.it/micciolo/PeM/01-rng.pdf
# -------------------------------------------------------------------
# se si vogliono riprodurre i risultati dei Laboratori
# con le versioni pił recenti di R va eseguita la seguente istruzione
# -------------------------------------------------------------------
RNGkind(sample.kind="Rounding")
# -------------------------------------------------------------------


hat <- c(1:60)

set.seed(654321)
sample(hat, size=2, replace=FALSE)


nrep <- 100000
out <- matrix(nrow=nrep,ncol=2)
set.seed(123456)
for (i in 1:nrep) out[i,] <- sample(hat, size=2, replace=FALSE)

job <- function(x,i) sum(x == i)

ris <- apply(out,1,job,1)
table(ris)

ris <- apply(out,1,job,50)
table(ris)




# lasciato al lettore (pag. 583)
y <- numeric(60)
for (i in 1:60)
{
    ris <- apply(out,1,job,i)
    y[i] <- table(ris)[2]
}
y <- y/nrep

sort(y)
mean(y)
sum(y)

tmp <- cbind(c(1:60),y)
o <- order(tmp[,2]); tmp[o,]